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Abstract 

Explicit formulas for the Gauss decomposition of elliptic Cauchy type matrices are 
derived in a very simple way. The elliptic Cauchy identity is an immediate corollary. 



1 Introduction 



The matrix 



1 ^ 



has determinant 



C=[— , q ,reC\ (1) 

, c , = Ul<i<j<N(Qi ~ Qj)( r 3 ~ r i) ^ 
ril<jj<Af(% ~ r j) 

This identity was first obtained by Cauchy. Now well-known as Cauchy's identity, it has found 
applications in harmonic analysis, soliton theory, and relativistic Calogero-Moser systems. 

Its elliptic generalization involving Weierstrass' sigma function a(z) is less widely known. 
It is given by 

det f a(q i -r j + X) \ N = <r(A + Ek=Mk - r k )) IIi<i<j<jv ~ <lj)°( r i ~ r j) (3) 

This elliptic Cauchy identity dates back to a paper by Frobenius j3]. Like (j2J), it has shown up 
in various contexts, giving rise to different proofs, cf . Refs. [TJ [6j [5j H] . 

Clearly, the identity applies to any minor as well. Moreover, after multiplication from the 
left and right by diagonal matrices (leading to so-called Cauchy-like matrices) one can still 
evaluate minors explicitly. 

Our perspective, which eventually led to this note, stems from the study of the Lax matrices 
of Calogero-Moser type systems. In particular, we wished to find the Gauss decomposition of 
the elliptic Cauchy-like matrix CV(A) given by (J5J) below, i. e., to represent it as 

Cjv(A) = UDL, (4) 

where U, D and L are upper-triangular, diagonal and lower-triangular matrices, respectively. 

In principle, this decomposition can be obtained by invoking two previously known results. 
Specifically, the Frobenius formula (J3J) can be combined with a theorem saying that the elements 
of the relevant upper- and lower-triangular matrices can be expressed in terms of appropriate 
minors of the matrix to be decomposed [8]. Indeed, as already mentioned, the minors of an 
elliptic Cauchy-like matrix also follow from the Frobenius formula. 

In this note, we wish to report an alternative method to decompose Cjv(A) which we find 
interesting and insightful. First of all, it is very economic, inasmuch as an exposition of the proof 
of the general Gauss decomposition formula and whichever of the known proofs of the Frobenius 
formula would require far more space and time. Secondly, our direct Gauss decomposition leads 
to a remarkably simple new proof of the Frobenius formula itself, and also reproduces some 
other results of interest as easy consequences. 
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2 The decomposition formula 



Our proof of the following decomposition formula is self-contained, except for its use of the 
3-term identity of the cr-function, 

a(z+a)<j(z— a)a(b+c)a(b— c)+a(z+b)a(z— b)a(c+a)a(c— a)+a(z+c)a(z— c)a(a+b)a(a— b) = 0. 

(5) 

We recall that this identity follows directly from the well-known relation between the Weierstrass 
p-function and the cr-function, 

/ \ / x a(y + x)a(y - x) 

p{x) - p(y) = — rWfT2 — ' 6 

cf. [7]. (Indeed, one need only divide ([5]) by a(z) 2 a(a) 2 a(b) 2 a(c) 2 and use (Q.) 
Theorem. Let qi, q^, n, r^, A be complex variables and introduce 

N 

\ N = A, A fe _i := \ k + q k - r k ee A + ^(<?j -Tj), k = 1, . . . , N. (7) 

j=k 

Define the elliptic Cauchy-like matrix Cat (A) by 

<%M - ( n aJ F^) 1l7 i + x \ ( n aJ F^\ « = 1 

Vfcdii - Qk) J <r(X)a{qi - Tj) a(n - rj )J 

where is is understood that Ylk=N+i ■■■ = 1- ^ eri ^ e decomposition (j3J) given fey 

(£)«)-! = t/« = i« = Cf(Ai); U* = &>(\ j ), i<j; L ij = C^(Ai), i > j. (9) 
Proof. Substituting 

2z = qi-rj+q N -r N ,2a = qi+rj-q N -r N ,2b = qi-rj + q N -r N +2X,2c = qi~rj-q N +r N (10) 
in ()5]), we obtain an identity from which the first step 

n (\ \ - ( i n-i c N (\ N ) A / Cjv_i(Ajv_i) W l N _ x o A , , 

of an inductive decomposition follows by a straightforward computation. Here, Ijv-i stands for 
the unit (JV — 1) x (N — 1) matrix, c^ is a column vector whose (N — 1) components are Off, 
i = 1, N — 1, and 7^ is a row vector whose (N — 1) components are , j = 1, iV — 1. 

Applying the decomposition ffTTl) to the Cauchy matrix Cjv-i(A./v-i) and then to the Cauchy 
matrix Cjv_2(Ajv_ 2 ) etc., we arrive directly at the formula (TjJ with factors (jUj). □ 



AT, 
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It remains to discuss some consequences. First of all, the following result follows effortlessly. 
Corollary (elliptic Cauchy identity). The formulas fl4]) and (Q imply 

det(C„(A)) = n Ojf (A.) = ft a fH i = g( ,» + ^-'' ft " rt) ' - < 12 > 

which is just the identity on account of (jSj) . 

Secondly, it is worth noting that the decomposition formula provided by the Theorem 
remains valid if we replace the a-function by any non-zero odd holomorphic function that 
satisfies the 3-term identity fl5J). (Indeed, our proof only uses these properties of o~(z).) It is 
known [7] that all such functions are of the form 

a(z) = e a+ ^a(z), (13) 

where a and (3 are arbitrary complex numbers, and where it is understood that the rational 
and trigonometric/hyperbolic degenerations of the a-function are included. If we replace a(z) 
by the rational degeneration furnished by a(z) = z and take A to infinity, then we obtain the 
Gauss decomposition of the original Cauchy matrix (TjQ) as well as the determinant formula (J2J) 
from our result. 

Thirdly, from the trigonometric specialisation we can recover the Gauss decomposition of 
the Lax matrix of the relativistic trigonometric Calogero-Moser system that recently cropped 
up in the paper J2J written by two of us. In fact, it was our discussion of the decomposition 
of the latter matrix that eventually led to the general decomposition encoded in the above 
Theorem. 

Finally, we point out that a similarity transformation with the reversal permutation matrix 
can be applied to (jH) to obtain a 'lower-diagonal- upper' version of the decomposition formula. 
After a relabeling 

px, . . . ,p N ->■ p N , . . . ,pi, p = q,r, (14) 
this decomposition has a well-defined limit for — > oo, by contrast to the one in the Theorem. 
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